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Role of Fermi surface anisotropy in the study of gap anisotropy using magnetic-field-angle 
dependence of thermal oscillations in A y Fe2- x Se2 superconductors 
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We present a numerical study of the field-angle resolved oscillations of the thermal conductivity and specific 
heat under rotated magnetic field in the A y Fe 2 -xSe2 [A=K,Rb,Cs,(Tl,K)] superconductors, using realistic two- 
band Fermi surface parameterization. Our key finding is that even for isotropic pairing on an anisotropic Fermi 
surface, the thermodynamic quantities exhibit substantial oscillatory behavior in the superconducting state, even 
much below the upper critical field. Furthermore, in multiband systems the competition of anisotropics between 
two Fermi surfaces can cause a double sign reversal of oscillations as a function of temperature, irrespective of 
gap anisotropy. Our findings put severe constraints on simple interpretations of field-angle resolved measure- 
ments widely used to identify the angular structure of the superconducting gap. 

PACS numbers: 74.20.Rp,74.70.Xa,74.25.Uv,74.25.N- 



The identification of the symmetry of the superconducting 
(SC) order parameter is an important step toward unraveling 
the pairing mechanism in any novel superconductor. For iron 
pnictides, the presence of hole and electron pockets at the T 
and M points has led to the proposal of ,s ± -wave pairing (T]- 
4] due to interband nesting between them. However, the re- 
cent discovery of the layered high-temperature superconduc- 
tors A 2y Fe2-a?Se2, with A=K,Rb,Cs,(Tl,K), has challenged the 
consensus for the pairing symmetry and mechanism of super- 
conductivity in this class of materials 0. The iron-selenide 
family has a crystal structure similar to the iron-pnictide ma- 
terial BaFe2As2, but with hole pockets eliminated completely 
from the Fermi surface (FS) at the T point in the Brillouin 
zone, yet the SC transition temperature T c is comparable to 
that of iron pnictides. Various theoretical proposals have been 
put forward which support either the survival of s-wave pair- 
ing [6-8 ], the emergence of nodal d-wave gap [9, 10], or more 
popularly nodeless d-wave gap |6]|7l[TT][T2l. Indirect experi- 
mental evidence suggests isotropic pairing symmetry lfT5lfT7lL 
consistent with isotropic gaps reported in angle-resolved pho- 
toemission spectra QjOQjSl. Therefore, direct high-precision 
imaging of the structure of the gap function and the location 
of the nodes, if they exist, is required. An effective and ac- 
curate technique for measuring the angular structure of the 
bulk gap relies on probing thermodynamic properties in a 
rotating in-plane magnetic field. For cuprate, pnictide, and 
heavy-fermion superconductors, this technique has been used 
widely to identify the SC pairing symmetry by mapping the 
field-angle dependence of the thermal conductivity or specific 
heat onto the angular structure of the SC gap and its pairing 
symmetry |[Tll20lJ30l. 

In this Letter, we demonstrate that detailed knowledge of 
the FS topology and parameters is necessary for relating the 
nature of the oscillations to the nodes or minima of the gap 
structure. This is especially important for materials where the 
FS anisotropy is substantial, as is the case in layered iron se- 
lenides. To be quantitative and unambiguous about the shape 



of the SC gap, it is required to incorporate realistic FS topol- 
ogy, Fermi velocities, and density of states (DOS) at the Fermi 
level into self-consistent calculations of thermal properties. 

To accomplish this goal, we focus on layered iron-selenide 
superconductors and study the information embedded in the 
angle-resolved specific heat coefficient, 7 = C/T, and ther- 
mal conductivity, k, in a rotating in-plane magnetic field using 
realistic tight-binding dispersions derived from first-principles 
electronic structure calculations. The main results of our cal- 
culations are: (1) For purely isotropic pairing symmetry, mod- 
erate FS anisotropies of layered iron-selenide superconductors 
are sufficient to introduce field- angle-dependent oscillations 
in the specific heat and thermal conductivity over a significant 
range of temperatures and at intermediate to high magnetic 
fields in the SC state. We find an inversion of the oscillation 
pattern as a function of temperature, which shows that oscil- 
lations are not a simple consequence of the anisotropy of the 
upper critical field. Therefore not all such oscillations at in- 
termediate fields can be taken as proof of strong anisotropy 
in the SC gap. (2) For isotropic gaps on the FSs, the oscil- 
lations in 7 may change sign once or twice as a function of 
temperature. We identify the out-of-phase FS anisotropies be- 
tween bands as the source for two sign reversals. (3) Complex 
field-angle dependence of the specific heat and thermal con- 
ductivity for anisotropic FSs suggests that comparison of both 
quantities with material-specific theories is needed to identify 
the pairing symmetry and gap structure. 

Anisotropy in FS and SC gap. — In the iron selenides the Fe 
vacancy completely eliminates the hole pocket at the T point, 
and the FS consists of two concentric electron pockets at the 
M point in the 2-Fe unit cell. This picture follows from first- 
principles calculations [32] and photoemission spectroscopy 
[QlOIIl). Here, we use a first-principles derived tight-binding 
parameterization of the electronic dispersion [ 6 ] with a weak 
k z dispersion as input to obtain all necessary FS parameters 
for a self-consistent transport calculation. Cuts of the corre- 
sponding FSs are shown in Fig. 1(a) with calculated normal- 
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FIG. 1. (Color online) (a) Different colors show different electron- 
pockets at k z = within the tight-binding model of the 2-Fe unit 
cell |6|. (b) Polar plots of computed normal-state DOS at the Fermi 
level show the out-of-phase anisotropy between bands (same color 
as in (a)), (c) The tight-binding in-plane Fermi velocities at k z — 
vs. azimuthal FS angle cj) exhibit out-of-phase FS anisotropy in units 
of lattice parameter a. (d) Gap functions A n (</>) on FS1 (solid) and 
FS2 (dashed) for pairing symmetries considered. The nodeless states 
have negligible FS anisotropy. (e)-(f) Field-induced total SC DOS 
at T = vs. energy at two representative field angles a for s and 
d x 2_ y 2 pairing. Note the low-energy crossings in the DOS (arrows) 
related to the low-T sign reversals in the oscillations of 7 and k; we 
used H/H C 2 = 0.5 for s wave and 0.1 for d x 2_ y 2 wave. 



state DOS in 1(b) and moderately anisotropic Fermi velocities 
in 1(c). These figures demonstrate the out-of-phase in-plane 
anisotropies of the FS parameters on the electron pockets at 
the M points. 

We consider three nodeless gaps with s, s ± = cos k x a + 
cos k y a, and extended d xy = sin (k x a/2) sin (k y a/2) sym- 
metry, shown in Fig. 1(d). Since the FSs are centered around 
M=(7T, 7T, 0), and its equivalents, all three including d xy are 
nodeless on the FS |6|. As all nodeless gaps exhibit very 
similar behavior, and thus we show detailed results only for 
the isotropic pairings s. For the nodal SC gaps, we con- 
sider two pairings symmetries as d x 2_ y 2 = cos k x a — cos k y a 



and d xy = sin k x a sin k y a, with detailed results presented for 
d x 2_ y 2 pairing. The gap structure for each pairing on the FSs 
is demonstrated in details in the supplementary material (SM) 

on. 

Brandt -Pesch-Tewordt (BPT) approximation.— We solve 
the quasiclassical Eilenberger equation within the extended 
BPT approximation CHU [3 EH gOj to solve for the 
field-angle induced SC DOS, N n (u,kf;H), together with 
the self-consistency equations for the SC order parameter, 
A n (kf; H), and transport lifetime, r n (cj, fe/; H). Here H is 
the magnetic field applied at angle a with respect to the (100) 
direction, kf is Fermi momenta and n = 1 , 2 is the band in- 
dex. The transport lifetimes encodes the combined effects of 
impurity and vortex scattering. The BPT approximation im- 
plies that the DOS is obtained by averaging the normal quasi- 
particle Green's function over the unit cell of the Abrikosov 
vortex lattice. This produces quantitatively correct results near 
the upper critical field over the range 0.5H C 2 < H < H c2 
l39l l42l l43ll for isotropic gap, but it extends to low fields for 
nodal and strongly anisotropic gaps dEBU- The SC gaps 
are evaluated by solving the coupled BCS gap equations for 
A n (fc/; H) at each applied field H. We simplify the problem 
by considering interband pairing only and eliminate the pair- 
ing potential in favor of the bare transition temperature T c q 
within weak-coupling theory [2j. Details of the calculations 
are given in the SM l34l . 

Based on the above-mentioned self-consistent solutions, the 
specific heat, C, and thermal conductivity, k, are computed 
numerically from the solution of the Eilenberger equations (3j 
4]. However, to get a qualitative understanding of how the 
interplay of FS anisotropy and gap anisotropy contribute to 
the results of C = C\ + C2 and k = k\ + ^2, we write down 
the approximate low-T expressions: 






00 u 2 (N n (Lj,k f ;H)) FS 
^ 4T 2 cosh(cj/2T) 2 ' 



(1) 



duo 



uj 2 (v x n (k f ) 2 N n (u, k f - H)r n (uj, k f ; H)) FS 



2T 2 cosh(cj/2T) 2 



(2) 



The angle-dependent SC DOS is given by the unit- 
cell averaged quasiclassical retarded Green's function g n : 
N n (uj,k f ;H) = -N fn (kf)Im g n (u>,kf',H)/ir. Here 
(. . .}fs stands for the FS integrals and v n is the Fermi ve- 
locity in each band, see Fig. 1(c). In the normal state the 
DOS becomes N n (u,k f ;H) = N fn (k f ) - l/|v n (fe/)| and 
r n (cj, kf,H) = r irnp ESI. 

The dominant contribution to the anisotropy in C n at low T 
originates from the anisotropy in the SC DOS at uj = 0, 



(N n (0,k f ;H)) FS 




Nfn(kf) 



/ 2AA n (fe /; ff) V 
V \v£(k r ,H)\ J 



, (3) 



FS 



where A = (hc/2\e\H) ' is the magnetic length of order the 
coherence length £ between 0.bH C 2 < H < H C 2, \v^\ is the 



component of the rescaled Fermi velocity normal to H, and 
A n is the impurity renormalized order parameter 0. For a 
cylindrical^ symmetric FS, the angle-dependence of \v„\ is 
determined solely by the field direction [3 ], its interplay with 
the profile of A n (fc/; H) gives the anisotropy of C n (a). For 
complex FSs, there is an additional weighting of the integral 
due to momentum-dependence of Nf n (kf) and v n (kf), lead- 
ing to strong field-angle oscillations even for an isotropic gap. 

Calculations for single-band models El El El S3 and ex- 
periments on several classes of materials (29j 08] 09) demon- 
strated that the anisotropy in heat capacity undergoes inver- 
sion as T and H change. We qualititively reproduce the gen- 
eral sign reversal of the oscillations for nodal pairings even 
after replacing quasi-cylindrical FSs with more realistic and 
material-specific FS anisotropies. In addition, for an isotropic 
gap, we find that even for a single strongly anisotropic FS, 
one or more inversion(s) of the oscillation can occur. Multi- 
band effects add additional complexity due to competing FS 
anisotropies and self-consistently evaluated multiple gap am- 
plitudes, and the intuitive one-to-one mapping between oscil- 
lations and nodal directions becomes easily lost at finite tem- 
peratures and fields. 

Results and discussions.— Figs. 1(e)- 1(f) show the field- 
induced SC DOS as a function of quasiparticle energy below 
the SC gap for a = 0° and a = 45° for s and d x 2_ y 2 OPs. 
We immediately see that the difference between SC DOS at 
these two representative angles changes sign at finite uj for 
both cases, opening the possibility for the low-T sign reversal 
of the oscillations in the specific heat as a function of temper- 
ature. 

We present the full angle-dependent profiles of C (a) and 
K,(a) for several temperatures at a representative low field for 
an isotropic s-wave gap (at H/H C 2 = 0.5) and a nodal d x i_ y i 
gap (at H/H C 2 = 0.1) in Fig. 2. It is interesting to note 
that there is a prominent angle dependence for the s-wave gap 
(Figs. 2(a) and 2(b)), whose nature is far more complex than 
what can be interpreted by conventional harmonics of pair- 
ing symmetries. Especially, at low T the peak position of 
C(a) is shifted from high-symmetry values and lies some- 
where between a = 0° to 45°, see Fig. 2(bl). Such com- 
plex field-angle dependence is a manifestation of the out-of- 
phase anisotropies on both FSs, shown in Fig. 1(a)- 1(c). For 
nodal d x 2_ y 2— pairing, the behavior of oscillations of C(a) 
and k(ol) is similar to results obtained for quasi-cylindrical 
FSs 13, however the amplitude of oscillations and the loca- 
tion of sign reversals are modified (56). 

We estimate the amplitudes of the fourfold oscilla- 
tions by defining C 4a (T) = U§ - n£ 5 , where n£ = 
[C(a,T)/T\/[C N /T c ] and «^(T) = [Tiff + ng ]/2 - n&, 
where II£ = [K xx (a,T)/T\/[Kff/T c ], and C N and k n are 
their corresponding normal-state values at T c 1571 . Such def- 
inition removes any twofold contribution from k. The cor- 
responding results are plotted in the lower panels of Fig. 2. 
We obtain several sign reversals in C/± a (T) and ^4 a (T) for 
both isotropic s- wave and nodal d-wave gaps. Earlier such 
sign-reversal feature was only found for highly anisotropic or 
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FIG. 2. (Color online) (al) Normalized thermal conductivity k/T 
along (100) direction (normalized to its normal-state value kn/T c ) 
as a function of in-plane field- angle a at fixed H/H C 2 = 0.5 for 
s wave, plotted from low to high T (bottom to top curves). Each 
subsequent curve is shifted vertically by 0.01 for clarity. Each curve 
is colored by the amplitude of the fourfold oscillation given in panel 
(a2); a uniform color map is used for values below —0.01 and above 
0.01. (a2) The fourfold amplitude of k/T is plotted as a function of 
T. The vertical arrows in the bottom row depict the temperatures at 
which curves in top panels are shown. Panels (bl)-(b2): Normalized 
specific heat for same parameters as in (al)-(a2). Panels (a) and (b) 
are for purely isotropic s-wave case, while similar plots in (c) and (d) 
are shown for nodal d x 2_ y 2-wa\e symmetry at H/H C 2 =0.1. 



nodal gap structure CHS |29l 0S 09). For this realistic FS 
parameterization of the double layered iron selenide, we find 
indications of two sign reversals even for the s-wave gap. Al- 
though the second sign reversal at high T may be difficult to 
discern, it is visible in Figs. 2(bl) and 2(b2) and as a white 
region in Fig. 3(al). We verified that the magnitude of oscilla- 
tions depends on the out-of-plane electronic hopping, defined 
in SM[34|. In addition, our calculations show that the ampli- 
tude of fourfold oscillations for s-wave pairing is roughly half 
of that for nodal pairing. 

It is noteworthy that the T-dependence of the fourfold os- 
cillations in Figs. 2(b2) and 2(d2) reflects on the energy de- 
pendence of the SC DOS in Figs. 1(e)- 1(f) at the same value 
of H. For example, for s-wave gap C(a = 45°) > C(a = 
0°) at low T in Fig. 2(b2), which corresponds to TV (a = 
45°) > N(a = 0°) at low energy in Fig. 1(e). The oppo- 
site anisotropy in C(a) at high T corresponds to the inver- 
sion of the DOS anisotropy at higher energies. We conclude 
that the anisotropy in C(a) for isotropic s-wave pairing is not 
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FIG. 3. (Color online) Contour maps of fourfold amplitude oscil- 
lations of normalized specific heat, C^a (top row), and normalized 
thermal conductivity K^a (bottom row) in the H-T phase diagram. 
Each column denotes a different gap symmetry. All plots use the 
same color map (red to blue); a uniform color map is used for values 
below -0.01 and above 0.01. Note the fourfold amplitude is given 
with respect to H//(100), i.e., a negative value corresponds to a 
minimum at a = 0°. Here T C (H) is defined by the vanishing of 
both gaps for each symmetry, which determines the line of the upper 
critical field H C 2- 



merely a manifestation of the anisotropy in H c2 (irrelevant at 
low fields, since it is tied to the FS shape), but is a reflection 
of the field-induced spectral-weight redistribution inside the 
gap. 

In Fig. 3, we show the contour map of the amplitude of the 
fourfold oscillations extracted from the normalized 7 = C/T 
(top row) and k/T (bottom row) for two nodeless and two 
nodal gaps. Earlier calculation using quasi-cylindrical FSs 
showed that the specific heat oscillation simply changes sign 
between the d xy and d x i_ y i symmetries, while the overall 
phase diagram remains very much the same between theml3l. 
After the inclusion of realistic and material- specific FSs in 
this work, we find substantial quantitative differences in the 
location of the sign-reversal lines between these nodal gaps 
in Figs. 3(cl) and 3(dl), due to the interplay of the SC order 
parameter with the FS anisotropies. Hitherto unknown is the 
intriguing result of both low-T (strong) and high-T (weak) 
sign reversals in the fourfold oscillations of C (a) and k(ol) 
for isotropic gaps, at moderate and high magnetic fields, see 
Figs. 3(a) and 3(b). We verified for s-wave pairing that the 
high-T sign reversal is robust and remains at nearly the same 
location for a single-band superconductor with identical FS, 
while the low-T feature disappears. Similarly, the low-T sign 



change does not exist for two-band models with similar (in- 
phase) angular variations of Fermi velocities. 

The striking feature of the phase diagrams for the heat ca- 
pacity for nodal and nodeless cases in the top panels of Fig. 3 
is that they all look qualitatively similar in the sense that they 
all exhibit sign reversals around the same T and H. The same 
is true for the phase diagram of the K,4 a term in the thermal 
conductivity. Furthermore, the sign of the fourfold oscilla- 
tions can sometimes be different for C± a and K^ a in the same 
region of phase diagram, as in the high-H and low-T region 
and vice versa for d x i_ y 2 -pairing in Fig. 3(d). This suggests 
that field-angle studies of each quantity alone are insufficient 
to distinguish between pairing symmetries. A simultaneous 
study of both C(a) and K,(a), including the comparison of 
the complex angle-dependent profiles and T behavior, is nec- 
essary to image the gap structure. Our field- angle-dependent 
results of the nodal d x 2_ y 2-wave gap in the H-T phase di- 
agram are in qualitative agreement with recent specific heat 
dataofCeMn 5 |49l . 

Conclusions. — The main conclusion of our work is that a 
mere observation of oscillations and sign reversals in C(a) 
or ft (a), combined with the H c2 anisotropy, is insufficient 
to identify the presence of nodes or minima in the gap, and 
their interpretations require detailed knowledge of the un- 
derlying FS anisotropy. For multiband systems, the situa- 
tion is further complicated by the interplay between multi- 
band FS anisotropies and multiple SC gaps in that substan- 
tial fourfold oscillations and sign reversals can occur even for 
purely isotropic gap. These results are robust in the region 
of 0.bH C 2 < H < H c2 for isotropic gaps, and its region of 
validity extends to lower field with increasing gap anisotropy. 
Our results suggest that not only realistic theoretical calcula- 
tions including field-induced impurity effect on mutiband sys- 
tems are necessary, (53) the field-angle measurements should 
also be compared with T-dependence of the penetration depth, 
specific heat, and residual electronic term of n/T measure- 
ments for the detection of pairing symmetry! 50-52]. In fact, 
in other probes such as quasiparticle interference (QPI) pat- 
tern seen in scanning tunneling microscopy/ spectroscopy 
(STM/S),|55| the consideration of the field induced impurity 
effect and the inclusion of realistic FS anisotropy should play 
an equally important role in the interpretation of data|54|. 
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SUPPLEMENTARY MATERIAL 

A detailed account of the extended Brandt-Pesch-Tewordt 
(BPT) quasiclassical approximation in the vortex phase has 
been given in a series of papers |1-5|. In this supplementary 
material we summarize the salient ingredients of this approx- 
imation and the numerical details of our quantitative studies. 

Fermi surface parameterization.— In order to calculate 
thermal properties in rotating fields, we use as input Fermi sur- 
faces (FSs) and dispersions of A y ¥Q2- x ^2 derived from first- 
principles electronic structure calculations to obtain an effec- 
tive tight-binding model in the 2-Fe unit-cell notation given 
in Refs. [6.7]. In accord with experimental observations, and 
essential for self-consistently calculated Green's functions for 
in-plane field rotation, a weak k z dispersion is added with 
hopping parameter t z ~ 0.3£, where t is the in-plane nearest 
neighbor hopping. Finally, a FS parameterization is obtained 
from the effect tight-binding model. For each FS integration, 
we used 800 points along the Fermi line in the {k xi k y ) plane 
and 9 k z slices in the Brillouin zone. We confirmed numer- 
ical convergence by checking that the results do not change 
when we increase the number of k points; this test was done 
for several points in the T-H plane. 

Multiband gap symmetries.— The superconducting (SC) 
pairing functions of all order parameters (OPs) considered for 
AyFe2-xSt2 are given in the main text. Fig. SI shows the 
two-dimensional view of of the corresponding OP on each 
Fermi surface sheet. The conventional s-wave gap is isotropic 
at all momenta as shown in Fig. Sl(a). The extended d xy 
order parameter is nearly isotropic on each of the FSs, but 
changes sign between FS sheets located at different M and 
equivalent points in the Brillouin zone (BZ), see Fig. Sl(b). 
This is one of the candidate order parameters for this system, 
which gives rise to a spin resonance (51 0. The s ± gap has a 
nodal line at the mid-point between the T to M direction, but 
that nodal line does not touch either FS. Importantly, due to 
the absence of a hole pocket at the T point, there is no sign 
change in the SC gap in A v Fq2-x^2 for this pairing. The 
nodal line passes through the zone boundary and zone diago- 
nal for d xy and d x 2_ y i -pairings as presented by dashed lines 
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FIG. 4. FIG. SI. (Color online) Order parameters and phase diagrams for the model pairing states in A y Fe2-xSe2. Panels (a)-(e): FS maps of 
the order parameter values for various pairing symmetries as indicated in each Panel. As shown in the color scale on the left, red means negative 
order parameter (OP), white denotes nodal regions, and blue indicates positive OP. The dotted lines depict the nodal line. Panels (f)-(g): Phase 
diagram of the fourfold oscillations for s ± -pairing extracted from heat-capacity, C^a, and thermal conductivity, K^a, respectively. 



in Figs. Sl(d)andSl(e). 

Computational details.— For each pairing symmetry, the 
coupled order parameters are computed self-consistently at 
each given magnetic field (H applied at angle a to the (100) 
direction, and temperature (T). To plot the phase diagram 
in Fig. 3 in the main text, we took a mesh of 35 field points 
between zero and H C 2, 100 temperature points from zero to 
T c0 , and used 31 a points from zero to 90° to extract the 
anisotropic terms in the heat capacity and the thermal con- 
ductivity. The phase diagram for s ± pairing is very similar to 
that of the nodeless s and d xy gaps, see Fig. Sl(f)-Sl(g), and 
thus not included in the main text. For all calculations, we 
considered purely intraband impurity scattering in the clean 
limit, l/r irnp = 0.01 x 27rT c0 , where T c0 is the bare transi- 
tion temperature and the scattering phase shift is chosen to be 
S = 7r/2 (unitarity limit). 

Methods and formulas.— If the Fermi velocity in band n 
is denoted by v n (kf), the corresponding normal-state DOS at 
the Fermi level is Nf n (kf) ~ l/|v n (fe/)|. The wavevector 



kf lies on the respective FS. When a magnetic field H is ap- 
plied along an angle a with respect to the x-axis, the relevant 
parameter in the SC state is the component of the Fermi ve- 
locity normal to the applied field, which, in energy units, be- 
comes v fn ((j), H) = v^(k f )/2A, where A = (Hc/2\e\H)^ 2 
is the magnetic length of order of the intervortex distance, <j) 
is the FS angle with respect to the k x axis, and v^(kf) is the 
rescaled component of the Fermi velocity perpendicular to H. 
The unit-cell averaged Green's function at the fermionic Mat- 
subara frequency, uj v , follows the notation of Eqs. (46)-(48) 
inRef.El 



g n (iuj u ,k f ;H) 



(4) 



I ioj u - 



A n (k f ;H) 



-Yl n (iuj u ,kf,H) 



Here A n and uj u are the order parameter and Matsubara fre- 
quency renormalized by the impurity self-energy, which is 
evaluated in the T-matrix approximation. The self-energy H n 
is given by (iuj u - D n )~ 2 = i^/7rW^(iuj/v fn ). W'(z) is 



the first derivative of the complex- valued function W(z) = 
exp (— z 2 )erfc(— iz). In contrast to the Doppler shift approxi- 
mation, both the real and the imaginary parts of S n contribute 
to the SC DOS, and their interplay as a function of energy, H 
and T, determine the sign reversal in the fourfold oscillation 
of the SC DOS. These effects have been extensively studied 
earlier using a single quasi-cylindrical FS and nodal gap, and 
a minimal 2D model for two-band systems, see for example 
Refs.EQH 

The order parameters are calculated self-consistently from 
the coupled gap equations of the two-band model 



A n (k f ;H) = 



(fc/, k s ')N fn , (k/)f n , (iu y , fc/; H) ) (5) 



Here f n t is the anomalous Gorkov function (off-diagonal 
Green's function). We used a factorized pairing potential at 
the Fermi surface as V nn '{kf,k/) = V nn >y n (<l))y n >(<l)'), 
with y n ((f)) the azimuthal angle dependence of the SC gap, 
see Fig. 1 in main text. For simplicity, we consider purely in- 

I 



terband pairing V12 = V21 = — V and eliminate V in favor of 
the bare transition temperature T c0 using weak-coupling the- 
ory 12. 

The specific heat, C = C\ + C2, and thermal conductivity, 
k = K\ + K2, are calculated from the approximate expressions 
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(7) 



We note that close to the transition the full expression for the 
heat capacity computed from the entropy includes the temper- 
ature derivative of the gap functions A n . However, inclusion 
of these terms only minimally affects quantitative aspects of 
the results away from the transition. The field-induced SC 
DOS in each band, N n /N n f = — Im g n /7r, is calculated us- 
ing analytical continuation iuj v — )■ uj+iS to obtain the retarded 
Green's functions, and the transport lifetime due to both im- 
purity and vortex scattering GlIHEJ: 
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(8) 



Since the function x 2 /cosh(x/2) 2 has a peak at x ~ 2.5T, 
the heat capacity at low temperatures predominantly probes 
the anisotropy in the total SC DOS, N(uj = 2.5T,fc/;iJ). 
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Using the expansion of the error function, we obtain two lim- 
iting values for W'(z): W'(0) = 2i/^ and W'(z > 1) w 
—i/^/irz 2 . Thus the SC DOS for each band n becomes 



N n (uj ] H) = {N n (uj,k f] H)) F s^{ 
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(9) 



where A n and S are the impurity renormalized order parame- 
ter and quasiparticle energy, respectively. The first line above 
only makes physical sense when the BPT approximation is 
valid at low energies, i.e., for nodal gaps. In that case at low T 
(low energy) and at low fields, where A n (fe/; H) only weakly 
depends on the direction of the field, the SC DOS depends 
predominantly on the orientation of \vf n (kf, H)\ relative to 
the minima of A n (fc/; H). At uj = the inversion of the SC 
DOS as a function of the field for nodal gaps can be obtained 
in analogy to Refs. 0E). 

At higher energies, the second line of Eq. ^ gives the BCS 
result for the SC DOS and therefore field-angle variation en- 
ters via the anisotropy of the upper critical field that influences 



A n (kf ; H) in the vicinity of the transition. This result is valid 
for both nodal and nodeless gaps, including the fully isotropic 
situation. Crucially, for anisotropic Fermi surfaces the H C 2 
anisotropy in the order parameter is weighted by the normal- 
state DOS, Nf n (kf), leading to a complex behavior including 
the switching of the minima and maxima found in our Let- 
ter. In this regime, however, the energy width of the Fermi 
weighting factor in the integral exceeds the gap amplitude and 
a full numerical evaluation is required. The results of such a 
self-consistent analysis are presented in the main text. All our 
results are consistent with the general observations based on 
this expansion. 
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